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ABSTRACT 

A real time, man controlled simulator exhibiting the 
characteristics of a 3000 Ton Surface Effect Ship ( 3K-SES ) 
was developed using a set of simplified 5 degree of freedom 
equations of motion. 

A software program was developed to compute and present 
to the SES pilot on a real time graphics output "alpha- 
numeric-perspective" visual cues pertaining to the vehicle 
states . 

Man generated control responses during 3K-SES accidental 
effector nozzle failure tests and the requirement for 
computer assist control systems during the recovery phase of 
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NOMENCLATURE 



A wp 

A 31 

A 33 

A w2 

A wl 

A 22 

6 

C DX 

C DY 

C DZP 




2 

Average sidewall wetted area, ft 

starboard side 

. 2 
Average sidewall wetted area, port side ft 



Added mass coefficient in roll force ft-slug 

equation 



Added mass coefficient in pitch force ft-slug 

equation 



Average wetted sidewall area of the bow 

Average wetted sidewall area of the stern 

Added mass coefficient in yaw force 
equation 

Sideslip angle 

Coefficient drag on x-direction 

Coefficient drag in y-direction 

Sidewall roll moment lumped parameter 
coefficient 

Bow pitch force lumped parameter 
coefficient 

Sidewall roll force 

Sidewall starboard buoyancy force 

Sidewall port buoyancy force 

Stern buoyancy force 

Bow seal pitching force 

Bow buoyancy force 

Gravitational acceleration 

Moment of inertia about x-axis 

Moment of inertia about y-axis 



ft' 



ft 



slug s lb-/ 
ft 



rad 

lb-s 2 lb /ft 2 
r m 

lb- s 2 lb /ft 2 
f m 

2 2 
lb f sVft 



non- 

dimensional 

lbs 

lbs 

lbs 

lbs 

lbs 

lbs 



ft/s* 
slug ft 
slug ft 



2 

2 
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p 



p' 

* 

p 

q 

q 

r 

r 



Moment of inertia about z-axis 
Track deviation performance index 
Control effort performance index 
Overall performance index 

Summation of moments about x-axis 

Length of sidewall 

Actual draft of port sidewall 

x-direction displacement of hull 
drag centroid 

Average draft of SES sidewall 

Average draft of bow seal 

Length from center of gravity to stern 

Average draft of sidewall 

Pitch moment lever arm for bow sidewall 
buoyancy force 

Pitch moment lever arm for bow seal force 

Mass of the rigid ship 

Summation of moments about y-axis 

Summation of moments about z-axis 

Plenum pressure 

Lumped drag centroid point 

Lumped drag centroid point 

Roll acceleration 
Pitch rate 
Pitch acceleration 
Yaw rate 

Yaw acceleration 



slug ft 2 
-. 2 . 

ft sec 
,2 

rad sec 
non- 

dimensional 

lbs-ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

lbs-ft 
lbs-ft 
lbs ft 2 
non- 

dimensional 

non- 

dimensional 

rad/sec 2 

rad/ sec 

rad/sec 2 

rad/sec 

rad/sec 2 
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S 1 


Turning moment lever arm of no. 1 engine 


ft 


S 2 


Turning moment lever arm of no. 2 engine 


ft 


S 3 


Turning moment lever arm of no. 3 engine 


ft 


s 4 


Turning moment lever arm of no. 4 engine 


ft 


T 7 


Total thrust magnitude on no. 1 engine 


lbs 


T 8 


Total thrust magnitude on no. 2 engine 


lbs 


T 9 


Total thrust magnitude on no. 3 engine 


lbs 


T 10 


Total thrust mangitude on no. 4 engine 


lbs 


T 

forw 


Total forward thrust vector of effectors 


lbs 


T . 
side 


Total side thrust vector of effectors 


lbs 


T 

yaw 


Total turning moment generated by 
effectors 


lbs ft 


u 


Velocity in x-direction (surge) 


ft/sec 


• 

u 


Acceleration in x-direction 


ft/sec 


V 


Velocity in y-direction (sway) 


ft/sec 


• 

V 


Acceleration in y-direction 


ft/sec 


V 

S . 


Total velocity 


ft/sec 


w e 


Width of bow seal 


ft 


X 


Summation of forces in x-direction 


lbs 


X o 


X coordinate of SES 

NAV 


ft 


• 

X 

o 


X velocity of SES 

NAV J 


ft 


Y 


Summation of forces in y-direction 


lbs 


Y o 


Y. TWT7 coordinates of SES 
NAV 


ft 


• 

Y o 


Y NAV velocit Y of SES 


ft 


<5 


Effector angle commanded 


rad 


6 7 


Effector angle of no. 1 nozzle 


rad 


6 8 


Effector angle of no. 2 nozzle 


rad 
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6 9 


Effector angle of no. 


3 nozzle 


rad 


6 10 


Effector angle of no. 


4 nozzle 


rad 


'P 


Heading angle of SES 




rad 


<t> 


Roll angle of SES 




rad 


e 


Pitch angle of SES 




rad 


p 


Density of sea water 




lbm/f t 
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I. 



INTRODUCTION 



Extensive analysis of the dynamics of the Captured Air 
Bubble (CAB) Surface Effect Ship (SES) has been accomplished 
in the past. Such analysis has produced equations of motion 
of the SES type vehicle in six degrees of freedom [Ref. 1] . 
Analysis at the Naval Postgraduate School has centered 
around the XR-3 and 100-B SES [Refs. 2,3], with current 
analytical effort being directed toward the 3K-SES [Ref. 4], 
The primary tool for study of the 3K-SES characteristics 
is the highly comprehensive five degree of freedom (DOF) 

SES DBS Program developed by the David W. Taylor Naval Ship 
Research and Development Center [Ref. 5] . A version of 
this program used at NPS is called DBSIM5D. This data base 
program allows very detailed study to be conducted on the 
3K-SES response to precomputed control inputs such as step 
or ramp thruster commands and their associated accidental 
failure possibilities. 

A natural extension in the study of control of the 
3K-SES is to evaluate the vehicle's response to an on-going 
"man generated" control input on a real time basis. The 
rationale for developing a "man-in-the-loop" SES simulator 
would be to 1) test the man in responding to accidental 
failure situations that may be experienced during opera- 
tional conditions, 2) train personnel in operator control 
techniques unique to SES type vehicles, 3) provide a means 



14 



to test hardware to be utilized in SES craft, and 4) provide 
a design tool for development of control schemes or system 
structural changes to improve performance. 

Accordingly a simplified, non-linear five degree of 
freedom Real Time Simulator (RTS5D) for the 3K-SES was 
developed such that a continuous, on-going, observable, 
real-time solution is interfaced with man-generated control 
inputs . 
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II. EQUATIONS OF MOTION 



A. COORDINATE SYSTEMS AND ASSUMPTIONS 

It was recognized that the RTS5D would require accurate 
and yet simplified equations of motion when compared to 
those of the more extensive DBSIM5D. The rationale for 
requiring simplification was to start with a right hand 
side (RHS) that could be utilized effectively in a real 
time simulation. To satisfy this requirement a model was 
selected that would use point source of drag forces con- 
sisting of only a few lumped-coefficient terms. The 
equations used are based upon those developed in Ref. 1 
which follow this concept. Additional assumptions used in 
the simplified 5 DOF equations were as follows: 

1) All accelerations are measured at the center of 
gravity. 

2) Vehicle is "f ree-to-heave" . 

3) All cross coupled moments of inertia are zero. 

4) All cross products of angular velocities in force 
equations are zero. 

5) All roll and pitch angles used in force equations 
are subject to small angle approximations. 

6) Calm water conditions. 

7) Mass and mass distribution of vehicle are constant. 

8) Effect of changes of aerodynamic force are neglected. 
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Utilizing the above assumptions, the simplified 5 DOF 
equations of motion are developed using the coordinate 
systems shown in Fig. 1 and Fig. 2, the force vector 
diagrams shown in Fig. 3, Fig. 4A, Fig. 4B, and the 
following equations from Newton's laws of motion: 



SURGE 


m (u - 


vr ) 


= X 


SWAY 


m (v - 


ur ) 


= Y 


YAW 


x z * 




= N 


PITCH 


I y q 




= M 


ROLL 


I.. p 




= K 



X 



where : 
m 
u 
v 
r 



P 

q 



I 

I 

I 

X 

Y 

N 



Z 

Y 

X 



mass of the rigid ship 

velocity in x-direction (SURGE) 

velocity in y-direction (SWAY) 

angular velocity about the z-axis 

angular velocity about the x-axis 

angular velocity about the y-axis 

moment of inertia about the z-axis 

moment of inertia about the y-axis 

moment of inertia about the x-axis 

summation of forces in x-direction 

summation of forces in y-direction 

summation of moments about the 
z-axis 



, slugs 
, ft/sec 
, ft/sec 
, deg/sec 
, deg/sec 
, deg/sec 
, lb -ft 2 

' lb m - ft2 
' lb nf ft2 




, ft-lb f 
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M = summation of moments about the 

y-axis , ft-lb f 

K = summation of moments about the 

x-axis , ft-lbf 



Additionally the following navigation relationships 
are utilized (see Fig. 1) . 

X q = u cos ip - v sin ip 

Y q = u sin \p + v cos tp 

2 2 

V = u + v 

s 

6 = tan 1 (-v/u) 



where: ip = heading angle 

6 = drift angle 

6 = thrust vector angle 



Also, 



X = / X. dt 

o J o 

Y = f Y dt 

o ; o 

r = / r dt 

P = / P dt 

q = / q dt 



u = / u dt 

v = | v dt 

ip = / r dt 

<p = f p dt 

0 = / q dt 
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NOTE: 1) <J; = 0 when vehicle heads parallel to + 

axis. 

2) Vehicle in figure above shown with negative 
thrust vector angle 6, positive ip , negative 
3 , negative sway velocity v, i.e., in a 
right turn. 

Figure 1 

Definition Of Coordinate System (Part I) 
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NOTE: Direction for positive velocity and 

angular rates are shown 

Figure 2 

Definition Of Coordinate System (Park II) 
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NOTE 





: 1 ) 
2 ) 

3) 

4) 



p and p' are equivalent point force centroids. 

, 6 0 , 6 n , 6, n are thrust vector angles and 

/ o y j. u 



are not required to be equal. 

, Tg, T , T q are thrust magnitudes and are 
not required to be equal. 

Force directions shown are for a right turn 



Figure 3 

Surface Effect Ship (Top View) 




NOTE: 



1) Right Turn Forces Acting, i.e., ^ , 6g 

6^^ are negative 

2) T SIDE = V in yV ln VV in V T 10 sin 



Figure 4A 

Surface Effect Ship (/-'tern View) 
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NOTE: 1) T f0rw =T 7 co S 6 7+ T 8 cos« 8 *T 9 cos6 9 +T 10 oos5 10 

(see Fig. 3) 



Figure 4B 

Surface Effect Ship (Side View) 



B. FORCES AND MOMENTS 



1. Surge Forces (See Fig. 4B) 

The forward acceleration was determined to be 
generated by the summation of thrust vectored in. the forward 
direction T,. . The resultant acceleration of the SES was 



assumed to be opposed by a retarding force exhibiting a 
"velocity-squared" characteristic. Additionally, Newton's 
laws of motion specified a "v*r" product that contributed 
to the drag component of the equation when the vehicle was 
in a turn (centrifugal force reaction term) . The simplified 
equation of motion for the surge acceleration follows from 
a summation of forces in the u direction as depicted in 
Fig. 4B. 



2. Sway Forces (See Fig. 4A) 

The side acceleration of the vehicle was analyzed 
to be the result of the summation of the thrust vectored 
parallel to the vehicle's y-axis and a retardation "velocity 
squared" phenomena such as that produced by the cross flow 
drag term of the sidewall as described in Ref. 6. This 
retardation force is augmented by a "u*r" product as 
specified by Newton's Law. 



forw 



u 





+ vr 



v 




ur 
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2 

Note the requirement to force the C v /m term to 

maintain the sign of the sway velocity such that both left 

and right turns may be computed. 

3. Yaw Moments (See Fig. 3) 

The yaw acceleration was a turning moment summation 

which was defined to be influenced by three primary forces 

operating over moment arms. It is of interest to note that 

the third term is an "added mass" term which is considered 

to operate on the same pressure point as the C DY v|v| term. 

The term T is a summation of moments generated by the 
yaw 

four thrusters as shown in Fig. 3. 



T 

yaw 



= s 4 T 7 cos <$ 7 + s 3 T s cos ^8 ~ S 1 T 9 COS ^9 

S 2 T 10 COS ^10 ” 5 'o T 7 Sin 5 7 ” 5 'o T 8 sin ^ 8 
£ o T 9 Sln 6 9 “ £ o T 10 Sin ^10 

C DY * w V M /I z + V 1 ! + K 22 w,l y/' 1 z 



4. Pitch Moments (See Fig. 4B) 

The pitch acceleration was assumed to be the summa- 
tion of moments generated by forward thrust, T f orw ' the 
buoyancy of the sidewalls of the SES, F^ and F^, and a 
vertical force generated at the bow of the vehicle, F^* 

This vertical bow force was defined to be a lumped parameter 
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term which modeled the reaction force due to plenum pressure 
acting against the bow seal. 



*d = 


= average draft of bow seal 


A 

W 1 = 


“ * x * + (U x tan 9)/2)£ x 


ro 

ii 


= ~ (U x tan 0)/2)£ x 


F i = 


= A £ pg 

W 1 x 2 


F 2 = 


= a„ 2 ^2 pg 


F 3 = 


' C dpV'e U d - *31 tan 9) 


£ 

X 1 


= average draft of sidewall 


P b = 


= plenum pressure 


A 

W 1 = 


= average wetted sidewall area of the stern 


£ x s 
X 2 


= width of one sidewall 


w = 

e 


= width of bow seal 


A 

W 2 


= average wetted sidewall area of the bow 


A 33 


= added mass coefficient 


P 


= water density 


g 


=32.3 ft/s 2 



= l ever arm °f bow seal 
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lever arm of buoyancy force 



q = (T. H + F_J> + F»£ - 

^ forw p 3 31 23 



C DX U \ 



F l £ 3 ~ A 33 u< 3)/ i y 

Note the added mass term A^uq which was required 
to provide damping to the pitch moment and is specified 
in Ref. 1. 

It was found that the response of the SES to a 

pitch perturbation without the added mass component was 

undamped and approximately sinusoidal. Unlike the yaw 

equation where the added mass term was a small contributor 

to damping, the pitch added mass term was found to be of 

significant importance in modeling the known pitch motion. 

5. Roll Moments (See Fig. 4A) 

The roll acceleration equation was analyzed to be a 

summation of moments generated by buoyancy forces of the 

port and starboard sidewalls, F and F , the thrust 

sp ss 

vectored parallel to the y-axis of the vehicle, T , a side 

force C QY v|v|, and a lumped coefficient vertical force, 

F . The vertical force F was defined as the force 
sw sw 

generated in a turn due to the sidewall curvature (dead 
rise angle) acting against the cross flow of water. 
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sp 


pg 


A £ . 
wp dp 


ss 


pg 


A £ _ 

ws ds 





c 


__ V 1 V 1 


sw 


DZP ‘ ' 



= average draft of SES sidewall 



'dp 



'ds 



A 



wp 



A 



£. - £ tan d) 

d^ sw Y 



£. + £ tan <p 

d^ sw Y 



ws 



average wetted area port 
average wetted area starboard 



p=((F - F ) £ - T . , £ + C v I v I £ 

c sp ss sw side p DY 1 1 z 



- F £ -A... up)/I, 

sw sw 31 v i 



Again note the added mass term A^up/I^ which was 
found to be essential in modeling the damping phenomena of 
the roll motion. 

In summary, the simplified 5 DOF equations of motion 
used in RTS5D for the 3K TON SES are: 



SURGE u 



SWAY v 



T forw /m " C DX u2/m + Vr lFig ‘ 4BJ 
T side /m " C DY v l v l/ m ‘ ur [Fi 9- 4A J 
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YAW r 



PITCH q 



ROLL p 



T yai/ I z + C DY v l v l i t/ I Z + A 22 uv£ v/ 1 Z I Fl 9- 3] 



(T forw *p + F 3*31 + F 2*3 ‘ C DX u " 2 



~ ^1^3 ~ ^34 u< 5) /ly 



[Fig. 4BJ 



( (F - F ) £ - T . . £ + v I v I £ 

sp ss sw side p DY 1 1 z 



- F £ - A,- up ) /I 

sw sw 33 e 



X 



[Fig. 4A] 



D. PARAMETER IDENTIFICATION 

These equations are an extension of the 3 DOF flat turn 
SES model developed by Gerba and Thaler in Ref. 4. The 
identification of craft parameters C DX , C^ , and w is 
described in Ref. 4 and repeated here for completeness. 

The surge drag coefficient C is determined by selecting 
a steady state turn condition, where T f orw ' u > v » r > and m 
are known. 



u 



0 



T forw 

m 



C DX U 

m 



vr 



from which 



C 



DX 



T £ + 
forw 

2 

u 



mvr 
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The sway drag coefficient C Dy is determined by the same 
steady state turn condition which requires 



v 



0 



T 

sxde 

m 




ur 



where 



T . , - mur 

side 



'DY 



vM 



The sway drag moment arm follows utilizing the yaw 
acceleration equation 



r = 0 = 



Vw , C DY v l v l l w . A 22 uv 



which yields 



w 



-T 



'DY v|v| + A 22 uv 



The addition of the pitch and roll equations introduces 
additional lumped coefficients C^- and C . These are 

DP \J ur 

easily solved using known constants 1,1, C_„, Z , SL , 

J ^ w p DY z sw 

w , p, , and steady state values of F , F , F n , F 0 , 

w U U, b p b b 4 

T forw' T side' u ' v ' an< ^ r * icant to note that 

in a steady state condition the angular roll rate, p, and 
angular pitch rate, q, are both required to equal zero. 
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This is in contrast to the yaw rate, r, where a finite 
steady state value is desired in a turn. Thus the added 



mass terms of A 33 uq and A^up in the pitch and roll equations 
are not utilized in the determination of craft parameters; 
their function is strictly confined to damping of their 
respective accelerations. Therefore, utilizing the pitch 
equation with known steady state turn values. 



The roll equation under steady state turn conditions 
yields 



q 



0 




2, 31 tan 0) 




yields 



C 



DP 




P 



0 



( (F -F )l - T . , H + C-,, vlvl A 
' ' sp ss w side p DY z 




from which 



C 



DZP 
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where 



F = port sidewall buoyancy in pounds [Fig. 4A] 

sp 

F = starboard sidewall buoyancy in pounds [Fig. 4A] 

s s 

T . . = thrust vectored parallel to vehicles y-axis 

S1 e in pounds [Fig. 3] 

v = sway steady state velocity in feet/sec 

s s 

Cpy. = sway drag coefficient for steady state condition 

SL = sway drag moment arm in feet [Fig. 3] 

i = effector thrust moment arm in feet [Fig. 4A] 

P 

l = surge drag moment arm in feet [Fig. 4A] 

z 

= sidewall buoyancy force moment arm in feet 
[Fig. 4A] 
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III. INTEGRATION METHOD 



A. ERROR ANALYSIS 

A rectangular (RECT) type integration was used to 
compute the velocity and positional variables of the force 
equations. Selection of the type of integration over a 
more complex type such as Runga Kutta, Variable Step enabled 
the loop delay time due to integration to be held to a 
minimum . 

An analysis of the error introduced into the solution 
of the simplified 5 DOF equations through the use of 
rectangular integration vice Runga-Kutta, Variable Step was 
made. The simplified 5 DOF equations of motion were 
programmed on the NPS IBM 360/67 computer using the 
Continuous System Modeling Program (CSMP) . One series of 
computet runs was made utilizing rectangular integration 
over a spectrum of step sizes ranging from 0.01 secs to 
0,25 secs. A second computer run was made utilizing Runga- 
Kutta, Variable Step type integration. The steady state 
results are shown in Table I. Examination of the final 
value magnitudes established that no significant steady 
state error was due to rectangular integration for the time 
increments utilized. 

Transient error of the variables, r, <p , and 6 was 
examined. A series of computer runs using incrementally 
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larger rectangular integration step sizes was compared 
against an identical run using Runga-Kutta Variable Step 
type integration. The first peak absolute magnitude 
results are shown in Table II. The RTS5D had an iteration 
time of .084 seconds therefore the "transient" magnitudes 
were in error by approximately 2% to 3% due to rectangular 
integration. 
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DELT 

.050 

.075 

.100 

.125 



Note 



TABLE II. TRANSIENT FIRST PEAK ERROR ANALYSIS 
FOR RECTANGULAR INTEGRATION 



r 




0 


1.42% 


1.64% 


2.58% 


1.89% 


2.41% 


2.84% 


2.16% 


3.38% 


3.25% 


2.35% 


3.65% 


8.57% 



1) Step effector angle input of 25° 
and velocity (u) of 56 knts. 

2) All percentages are referenced to 
Runga-Kutta variable step magnitudes. 
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IV. IMPLEMENTATION 



A. PHILOSOPHY 

Given the decision to develop a real time motion analysis 
simulator for the 3K-SES, the various options available 
were examined. Prime requirements for such a simulator 
were the necessity for real clock interfacing, suitable 
input-output ports, and a dedicated Central Processing 
Unit (CPU) . 

The "all analog" computer simulation was particularly 
attractive from the aspect of interfacing man generated 
control inputs into an on-going real time solution. With 
the analog approach instantaneous solutions are obtained 
from patchboard amplifiers. Two considerations which did not 
lend themselves well to the all analog technique were how to 
maintain high accuracy solutions in wide ranging variables 
when constrained to a voltage limit of plus or minus 100 
volts, and how to output the system variables in such a 
manner that a man could physically interpret what was 
happening and respond with the required corrective, timely 
control effort. Additional supplementary problems associ- 
ated with the "quarter-square" multiplication error and 
amplifier drift/run duration limitations gave rise to 
examining other possible approaches. 

The "all digital" computer simulation offered high 
accuracy with wide ranging variables and a capability to 
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interface with two graphics display consoles. Additionally 
a line printer was available to provide "hard copy" of 
variables of interest as the solution was being processed. 
Two major disadvantages were the serial and discrete nature 
of the integration method of a digital computer, and the 
lack of high voltage, high accuracy potentiometer controls 
to allow the man to influence the motion of the vehicle as 
the solution was being presented to him. 

A desire to take advantage of the best features of both 
the analog and digital alternatives led to a decision to 
develop the real time simulator (RTS5D) of the 3K-SES 
utilizing the simplified 5 DOF equations of motion developed 
in Section II. 

B. REQUIREMENTS 

The following criteria were required of the RTS5D: 

1) The simplified equations of motion would be solved 
and the results output on a real time basis. 

2) The solution would be subject to real time control 
efforts generated by a man observing the output. 

Development of the simplified 5 DOF equations of motion 
and the use of rectangular type integration for simulating 
the 3K-SES ensured that loop delay attributable to digital 
computational needs would be minimal and thus the solution 
would be compatible to a real time environment. 

The requirement to output the results on a real time 
basis such that a man observed and responded with real time 



38 



control efforts necessitated the use of graphics terminals 
interfaced with a digital computer. Large volumes of data 
and graphical cues were presented to the pilot and updated 
on a real time basis in much the same manner as found on 
the bridge consoles of the 3K-SES shown in Ref. 7. Software 
algorithms designed to present "alpha-numeric" data and 
graphic perspectives corresponding to the vehicles current 
state variables enabled the observer to respond on a real 
time basis. 

Real time man generated control requirements necessitated 
the need for a hybrid environment and its associated inter- 
face between a digital computer and an analog computer. 

The RTS5D utilizes four (4) computer systems: One (1) 

medium size general purpose digital computer (XDS-9300) , 
two (2) small graphics digital computers (AGT-10) , and one 
(1) analog computer (CI5000) . The XDS-9300 has a CPU 
capable of handling 32K words of main (core) memory with a 
magnetic drum for secondary storage. The XDS-9300 may be 
programmed in assembly language or FORTRAN. The CI-5000 
is an analog computer with a full complement of digital 
logic. Hardwired trunk lines linking the XDS-9300 and 
CI-5000 allow hybrid computer operations. The two AGT-10' s 
are small general purpose digital computers with 8K of main 
memory and magnetic disk for secondary storage. The AGT-10 
systems have the capability of accepting input and of 
displaying output on their CRT screens. Both the CI-5000 
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and the AGT-10 have the capability of communicating with 
the XDS-9300 directly, but cannot communicate with each 
other except via the XDS-9300 [Ref. 8] . 

C. HARDWARE DESCRIPTION 

Hardware layout for the RTS5D is shown in Figure 5A. 

A "joystick" type high voltage potentiometer attached to 
the pilot's seat was hardwired to the CI-5000. This joy- 
stick served the purpose of the pilot's effector angle 
controller. Left or right stick transmitted an appropriate 
nozzle deflection angle in the form of a voltage level to 
the CI-5000 where it was quantized in an analog to digital 
process and introduced into the XDS-9300 's ongoing solution 
of the simplified 5 DOF equations. In like manner the 
voltage levels of six other "dial" potentiometers used to 
control thruster output levels and program modes were 
introduced into the solution process as deemed appropriate 
by the operator. The six potentiometers were located on a 
"thruster console" (See Fig. 5B) within reach of the pilot 
(See Fig. 5A) . Additional hardware included a standard 
digital computer line printer, an analog 8 channel strip 
chart recorder, and an analog X-Y plotter. 

D. SOFTWARE DESCRIPTION 

1. Introduction 

The RTS5D software package was a FORTRAN IV digital 
computer program which consisted of a main program and two 
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Figure 5A 

RTS5D Block Diagram 
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A - Thruster T-7 
B - Run, Hold, Restart 
C Thruster T-8 
D - Thruster T-9 

E - Accidental Failure Nozzle Angle 
Controller For Thruster T-10 

F - Thruster T-10 



Figure 5B 
Thruster Console 
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major subroutines. The program gave the operator the option 
to start, stop, and hold the dynamics solution process. 

The software package controlled a multiplexing algorithm 
used on the graphics display consoles and sampled all 
analog control and timing inputs. A program flow chart is 
shown in Fig. 6A and a detailed flow chart for the graphics 
multiplexing is shown in Fig. 6B. 

The program began by initializing all program 
variables and parameters and then requested from the operator 
any parameter changes. The request and input of changes 
were made on a TTY-35 Teletypewriter (See Fig. 5A) . The 
program then started the real time 1000 Hz clock timing 
circuits in the CI-5000 analog computer and entered the main 
program iteration loop. In the main loop the sampling of 
analog control inputs by the pilot was accomplished by 
subroutine ADDA. Craft force equations were computed and 
integrated. A perspective view of the simplified form of 
the SES and an electronic road were utilized in the RTS5D 
to give real time visual cues to the pilot on the status 
of the vehicle's states. This section of the display gave 
the pilot his primary guidance information and was updated 
on the pilot's graphics display console on every iteration 
of the main program loop. The program then evaluated 
whether an accidental failure mode had been selected by the 
operator. If the failure mode was selected the program 
searched the craft's states to see if the recovery criteria 
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Figure 6A 

RTS5D Program Flow Chart 
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DESCRIPTION OF INFO(I) 



INFO ( 1 ) u DYNAMICS 

INFO (2) u DYNAMICS 

INFO (3) v DYNAMICS 

INFO (4) ❖ DYNAMICS 

INFO (5) r DYNAMICS 

INFO (6) r DYNAMICS 

INFO (7) p DYNAMICS 

INFO (8) p DYNAMICS 

INFO (9) q DYNAMICS 

INFO (10) 4 DYNAMICS 

INFO (11) POSITION MAX VALUES 

INFO (12) POSITION MIN VALUES 



INFO (13) 
INFO (14) 
INFO (15) 
INFO (16) 
INFO (17) 
INFO (18) 
INFO (19) 
INFO (20) 
INFO (21) 
INFO (22) 
INFO(23) 
INF0(24) 



VELOCITY MAX VALUES 
VELOCITY MIN VALUES 
ACCELERATION MAX VALUES 
ACCELERATION MIN VALUES 
NAVIGATION DATA 
EFFECTOR/ THRUSTER DATA 
SEA TRACK NAV PLOT 
SEA TRACK NAV PLOT 
PERFORMANCE INDICES DATA 
PRESENT POSITIONAL VARIABLES 
PRESENT VELOCITY VARIABLES 
PRESENT ACCELERATION VAR. 



Figure 6B 

RTS5D Graphics Display Multiplex Algorithm 
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had been satisfied through the maneuvering inputs of the 
pilot. If the recovery criteria was fulfilled, the program 
went into a hold status and awaited further instructions 
by the operator on when to restart the run. If the recovery 
criteria had not been fulfilled the program continued the 
main loop by incrementing a graphics multiplex index and 
entering the graphics display multiplex algorithm to search 
for and display the appropriate information array associated 
with the current value of the graphics index I (See Fig. 6B) 
A graphics multiplexing algorithm was established to reduce 
the iteration time required for the main loop. Each informa 
tion array contained approximately 50 elements associated 
with positions on the display screens. Each array had a 
display time which was required to encode and project the 
array on the screen. This display time varied directly 
with the number of words in the array. The prioritization 
of the graphic display console information arrays was made 
such that the perspective of the SES and the electronic 
road were displayed on each iteration and were not subject 
to the program's graphics display multiplexing routine. 

All other information arrays were multiplexed on a time 
share basis such that on any main loop only two information 
arrays were projected on the screen, the SES/road and the 
array designated by the graphics display index I. This 
concept was mandatory due to the large volume of data/ 
graphics being presented. It is instructional to recall 
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the fact that all variables were computed and updated in 
the memory of the XDS-9300 on every loop, whereas the 
display routine was multiplexed on a time share basis. 

The main loop iteration time was 84 msec. Information 
arrays in the multiplex routine were updated on the graphics 
display consoles every 2.0 sec. On exit from the graphics 
display multiplex algorithm the program read the value of 
the real time clock and then reset the clock to zero. From 
this the A time magnitude was determined. Following this 
the program searches for print, restart, and hold flags 
which are set by the operator. The program loops back to 
the AD/DA subroutine and begins a new iteration. 

2 . Graphics Techniques 

The two AGT-10 Graphics Display Consoles were 
utilized to display real time data to the operator such that 
he could control the craft. The process of transmitting and 
displaying desired data to an AGT-10 was accomplished 
utilizing FORTRAN callable graphics subroutines in the 
XDS-9300. One graphics subroutine TEXTO was used to transmit 
alpha-numeric data and one graphics subroutine GRAPHO was 
used to transmit graphical data. 

Subroutine TEXTO transmitted alpha-numeric data to 
the AGT-10 screen. The positioning of the desired data on 
the screen itself was accomplished by utilizing a TEXTO 
screen coordinate system of dividing the 12 inch by 12 inch 
AGT-10 screen into 40 equally spaced lines from top to 
bottom with each line having room for 96 individual 
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characters on it as shown in Fig. 7A. Subroutine TEXTO 
statements analogous to FORTRAN "WRITE" and "FORMAT" state- 
ments were used to specify the letter character or variable 
magnitude to be displayed and their position on the display 
screen. 

Subroutine GRAPHO was used to transmit the "picture" 
portion of the output to the display screen. The coordinate 
system utilized by GRAPHO is a normalized X-Y system with 
the ordinate and abscissa intersecting in the center of the 
screen as shown in Fig. 7B. A point on the GRAPHO coordinate 
system corresponding to (l.,0.) would be located at the right 
edge of the screen centered between the top and bottom of 
the screen. A line could be generated on the screen by 
specifying the line start coordinates, line stop coordinates 
and commanding a "draw". The AGT-10 would connect the two 
designated points and a line would be displayed on the screen. 
Concurrently if only two point sources of light had been 
desired to be displayed, the command "jump" could be 
utilized such that the AGT-10 would not connect the second 
coordinate point with the first coordinate point. In this 
manner a series of commands was generated via subroutine 
GRAPHO whereby a picture was drawn on the AGT-10 screen 
using an array of sequential coordinate points and their 
respective "draw" or "jump" instruction. 

3. Analyst Graphic Display Console 

The display presented to the analyst during the 
conduct of a run is shown in Fig. 8. The magnitudes of the 
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VARIABLE STATUS 



POSITION 





X 


Y 


ZZ 


PSI 


PHI 


THETA 


Present 


2193.1 


-2172.0 


o 

• 

o 


1 

o 

• 

CT\ 


-.014 


. 005 


Maximum 


2193.1 


o 

• 

o 


o 

• 

o 


o 

• 

o 


.009 


.007 


Minimum 


o 

• 

o 


-2172.0 


o 

• 

o 


-0.6 


-.016 


.000 



VELOCITY 





U 


V 


W 


R 


P 


Q 


Present 


87.6 


4.86 


0.0 


-.01 


.002 


+0.00 


Maximum 


94.1 


5.18 


0.0 


0.0 


+ .161 


+ .110 


Minimum 


87.6 


0.0 


0.0 


-.012 


-.305 


-.034 



ACCELERATION 





UDOT 


VDOT 


WDOT 


RDOT 


PDOT 


QDOT 


Present 


-0.29 


-0.05 


0.0 


+ 0.00 


.001 


+ . 00 


Maximum 


0.0 


0.833 


0.0 


0.09 


.003 


+ . 00 


Minimum 


-0.36 


-0.06 


0.0 
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O 
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I- 1 


-.004 


-.00 



Note: Typical run values shown 



Figure 8 

Analyst Graphic Display Console 
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the variables were computed by the XDS-9300 on every iteration 
of the program. The XDS-9300 then followed a display multi- 
plexing algorithm (See Fig. 6B) which updated the variables 
every 2. seconds. The variables were categorized by position, 
velocity, and acceleration. Each category had three lines 
of data displayed. The "present" data line displayed the 
real time magnitude of the variable column it was located 
in. The "maximum" and "minimum" were the max and min 
magnitudes that had been computed for the variable during 
that particular run. At the start of a new run all variables 
are initialized (See Fig. 6A) . 

4 . Pilot Graphics Display Console 

The display presented to the pilot during the conduct 
of a run is shown in Fig. 9A and Fig. 9B. The pilot display 
contained three major sections, VARIABLE DYNAMICS, SEA TRACK/ 
EMERGENCY PERFORMANCE INDEX, CRAFT PERSPECTIVE/NORTH ROAD/ 

NAV DATA. 

a. Variable Dynamics 

The time histories of the accelerations and their 
respective velocities were presented in graphical form. An 
acceleration plot and its respective velocity plot shared 
the same electronic strip chart such that one was super- 
imposed over the other. The acceleration plots were 
displayed as "light" toned lines whereas the velocity plots 
were displayed as "heavy" toned lines. As the RTS5D program 
incremented in real time, it presented an incremental 
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Figure 9A 

Pilot Graphics Display Sections 
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Figure 9-B 

Pilot Graphic Display 
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increase in each displayed state time history. The time 
histories would progress in a lef t-to-right movement and 
upon reaching the right most limit of the VARIABLE DYNAMICS 
section were erased from the screen and restarted at the 
left side of the electronic strip chart. Execution of a 
maneuver of the vehicle would immediately be displayed in 
the time histories displays. 

b. Sea Track/Emergency Performance Index 

The sea track portion of this display area was 
a navigational plot of the ship trajectory (X o ,Y q ) (See 
Fig. 9B) . This track was updated every 2 seconds as imple- 
mented in the XDS-9300 multiplex algorithm. The track 
would erase itself and begin anew at the craft's current 
position every 3.5 minutes. 

The performance of the pilot during an acci- 
dental failure situation was measured and graded through 
the use of three performance indices (TRACK DEVIATION) , 

J 2 (CONTROL EFFORT) , and J 3 (OVERALL PERFORMANCE) . All 
performance tests were initialized such that the vehicle 
was on a heading of 000° (NORTH) and on the axis with 

a steady state velocity. 

J^, TRACK DEVIATION, was an "integral-error- 
squared-plus-time" performance index which penalized the 
distance that the pilot allowed the vehicle's c.g. to 
drift away from the axis in addition to the length 

of time the vehicle remained off the axis (See Fig. 10). 
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X NAV 




Heading=000 
Velocity= 
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Figure 10 

Track Deviation Performance IndexCJ^) 
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At the systems entry into an accidental failure mode of 
operation the XDS-9300 computed as 

*"final p 

J, = / (Y ) dt + TIME 

1 

J 2 / CONTROL EFFORT, was an index used to 
penalize the pilot for excessive control effort during the 
man controlled recovery maneuver segment of the accidental 
failure test. The J 2 index was an integral "error-squared" 
index, as 

t f inal 

J 2 = / (6 ) 2 dt 



Note that the 6 under the integral sign is the 6 associated 
with the nozzle deflection of the three operational thrusters 
of the SES, and does not include the nozzle angle of the 
failed thruster (i.e., starboard thruster). A graphical 
interpretation of J 2 is shown in Fig. 11. 

J^, OVERALL PERFORMANCE, was a magnitude 
representing the sum of J-^ and J 2 



57 



+5 (degrees ) 





Figure 11 

Control Effort Performance IndexC^) 
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At the completion of the test the performance algorithm 
in the XDS— 9300 was programmed to display the numeric 
scores representing and J 2 and the word evaluation, 
EXCELLENT, ABOVE-AVERAGE, or GOOD below the overall performance 
heading shown in Fig. 9B. The magnitude range of the 
performance index J 3 for each of the 3 word evaluations are 
shown in Table III. 

TABLE III. OVERALL PERFORMANCE 





INDEX 


RATINGS 




RATING 


EXCELLENT 


ABOVE AVERAGE 


GOOD 


RANGE 


J 3 < 250 


250 < J 3 < 500 


500 < J 3 



c. Craft Perspective/North Road/NAV Data 

Real time visual cues to the pilot were generated 
on the lower half of the pilot's graphics display console 
in the form of a perspective view of a simplified shape of 
the craft from the rear looking down and forward onto the 
craft and the area immediately ahead of the craft's path. 

An electronic north-heading road shown forward of the craft 
provided information to assist the operator in proper 
control of the vehicle. Numeric navigation data concerning 
the craft's position and the status of the nozzle angles 
and effector thrust settings were presented in the area 
immediately below the craft perspective (See Fig. 9A and 
Fig. 9B) . 
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The craft perspective view presented real time 
graphical information on the vehicle’s turning rate (r) , 
velocity (u) , pitch angle ( 0 ), and roll angle ( 4 >) . The 
electronic north road’s centerline was the positive X„, TI 
axis (See Fig. 1) . The north road was displayed on the 
pilot's screen whenever the pilot had maneuvered the craft 
to within 1500 feet of the X^^. axis and the heading of the 
vehicle was on a northerly course. At other times the north 
road was not visible to the pilot. When displayed, the 
north road provided the pilot with visual cues on the 
vehicle's heading and navigation position with respect to 
the fixed X NAV axis. When the pilot maneuvered the craft 
such that the north road was not visible graphical turn rate 
information was determined from a zero turn rate reference 
line located just above the bow of the SES. The pilot 
established a zero turn rate whenever the center point of 
the craft's bow was lined up just below the zero turn rate 
reference line. Variable speed lines moving from the 
craft's bow toward the stern gave a graphic representation 
of the forward velocity u (See Fig. 12) . A sequence of 
displays during a pilot commanded control effort of left 
nozzle deflection and then right nozzle deflection is shown 
in Fig. 13, Fig. 14, and Fig. 15. 

The north road positioning on the pilot's 
graphic display console was solved using the following 
assumptions : 
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CENTER OF ZERO TURN RATE 




NOTE: 1) Craft "shows center of bow is left of 

reference line hence a left turn status. 



Figure 12 

SES Perspective Information 
Without North Road 
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VARIABLE DYNAMICS SEA TRACK 





CONTROL INPUTS 




THRUST 


RUDDER 


DELT 


ELAPSED 


PITCH HEADING/RATE 


DRIFT 


ROLL 


SPEED 


400000 


.32 


.084 


10 


.52 


355/-1.0 


.02 


-.02 


56 


THRUST 1 


THRUST 2 


THRUST 3 


THRUST4 


RUD1 


RUD2 


RUD3 




RUD4 


100000 


100000 


100000 


100000 


.32 


.32 


.32 




.32 



NOTE: 1) Craft has been commanded by pilot into a left turn 
by deflecting nozzles .32 radians (+18°). Vehicle 
is entering left turn and is passing through a heading 
of 355° with a turn rate of -l°/sec (see NAVIGATION 
DATA readout) . 

2) SEA TRACK and electronic road show vehicle has dis- 
placed slightly left of north road centerline. 

3) VARIABLE DYNAMICS show velocity and acceleration time 
histories being generated as manuever develops . 

Figure 13 
SES MANUEVER (I) 
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VARIABLES 



SEA TRACK 





CONTROL INPUTS TIME 

THRUST RUDDER DELT ELAPSED 
400000 .32 .084 20.2 

THRUST1 THRUST2 THRUST 3 THRUST4 
100000 100000 100000 100000 



NAVIGATIONAL DATA 

PITCH HEADING/RATE DRIFT ROLL SPEED 
.58 320/1,2 .04 -.5 51 

RUD1 RUD2 RUD3 RUD4 

.32 .32 .32 .32 



NOTE: 1) Pilot has held left nozzle deflection angle of .32 
rad (+18°) and vehicle has displaced left of north 
road. Vehicle has heading of 320° or 40° off the 
north road centerline. 

2) SEA TRACK confirms trajectory of craft in left turn. 

3) VARIABLE DYNAMICS show +u,+v,-r, and q in steady 
state while roll rate p is still seeking steady 
state condition. 



Figure 14 
SES Manuever (II) 
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VARIABLE DYNAMICS 



SEA TRACK 





CONTROL INPUTS TIME 

THRUST RUDDER DELT ELAPSED 
400000 -.5 .084 33. 

THRUST 1 THRUST 2 THRUST 3 THRUST4 
100000 



NAVIGATION DATA 



PITCH HEADING/RATE 


DRIFT ROLL 


SPEED 


.51 


030/1.8 


-.3 .8 


50 


RUDl 


RUD2 


RUD3 


RUD4 


-.5 


-.5 


-.5 


-.5 



NOTE : 



1) Vehicle has been manuevered into right turn. 

2) Vehicle heading 030° on path which will intersect 
north road. 



Figure 15 

SES Manuever (III) 
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1 ) 



The view of the pilot was restricted to plus and 
minus 60 degrees of the heading of the craft. This 
corresponded to the far right and far left edges 
of the pilot's display screen, respectively. 

2) The starting edge of the road would remain fixed at 

a point ahead of the vehicle whenever the vehicle was 
positioned such that the pilot could see the road in 
his restricted viewing arc. 

3) The road would give visual cues on the craft's left- 
right displacement from the road centerline and on 
the relationship between the north heading of the 
road and the heading of the vehicle. 

4) The road's shape as it appeared to the pilot would be 

a wide path terminating at an apex point in the distance. 

A non linear relationship was used to compute the 
position of the north road (See Fig. 16A and Fig. 16B) . The 
Y NA v displacement of points A, B, and C on the graphics 
display console was made a function of Y q only. The Y NAV 
displacement of point D on the graphics display console was 
made a function of Y q and the arc tangent of the craft's 
heading. The X NAV displacement of the points A, B, C, and 
D were fixed at a constant value, thus recalling Fig. 7B, 
the computation of the coordinates of the points A, B, C, 
and D was 
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X NAV 




Figure 16 A 
North Road(I) 
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NAV 



Figure 16B 
North Road (II) 
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The displacement of points A, B, and C due to 
vehicle heading were neglected. 

E. OPERATIONAL MODES 

The RTS5D operated in any combination of modes as 
designated by the operator (See Fig. 17). Mode selection 
was accomplished via "on-off" switches located at the 
console of the XDS-9300 and at the thruster console located 
at the pilot's seat (See Fig. 5A and Fig. 5B) . The RTS5D 
ran in any combination of the modes available with the 
"RESTART" mode taking priority over all other modes selected. 

The FIXED THRUST and FIXED EFFECTOR ANGLE removed the 
appropriate man-generated control and utilized a constant 
value for T and 6 respectively. 

The PRINT mode allowed variables to be printed on the 
XDS-9300 line printer at intervals selected by the operator. 

The HOLD mode put the program into a loop whereby the 
display screens would continue to present the last image 
computed until such time as the mode was deselected. At 
this time the program would continue. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 



DESCRIPTION 



THRUST MAGNITUDE FIXED 
THRUST MAGNITUDE VARIABLE 
EFFECTOR DEFLECTION ANGLE FIXED 
EFFECTOR DEFLECTION ANGLE VARIABLE 
ACCIDENTAL FAILURE THRUST LOSS 
ACCIDENTAL FAILURE EFFECTOR REVERSAL 
PRINT 
HOLD 

RESTART ADJUST CONSTANTS IN PROGRAM 
RESTART RETAIN CONSTANTS 



Figure. 17 

RTS5D Modes of Operation 



69 



The RESTART mode re-initialized the program and then 
provided an opportunity to the operator to change the 
magnitude of any constant in the program. This feature 
allows the operator to carry out program modeling adjust- 
ments without the delays associated with re-compiling a 
modified digital program. Adjustments to constant magni- 
tudes are made via a peripheral teletypewriter (TTY-35) 

(See Fig. 5A) . 

The ACCIDENTAL FAILURE FFFECTOR REVERSAL mode deflected 
the nozzle angle of the starboard thruster (See Fig. 3) 

to 180°. This simulated the accidental reversal of an 
outboard thruster at full power. Simultaneously the 
program computed and displayed the performance indices 
discussed in the Software section of this report. Exit 
from this mode was accomplished when the pilot had maneuvered 
the craft back on track. At this time the program would 
"hold" while displaying the last computed values of all 
variables. The hold condition would continue until a 
restart was accomplished by the operator. 

The ACCIDENTAL FAILURE THRUST LOSS mode was entered by 
decreasing the thruster setting of any one, or combination 
of, thrusters (See Fig. 5B) . Exit from this mode was 
identical to the accidental failure effector reversal 
described previously. 

The FIXED STEP SIZE INTEGRATION mode removes the real 
time clock input and utilizes an operator specified time 
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increment on each iteration of the program loop. This is 
a non real time mode. The output displays and other modes 
are functional in this mode. 
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V. RTS5D RESPONSE CHARACTERISTICS 



A. VALIDATION 

The DBSIM5D program was utilized as the benchmark to 
validate the RTS5D simulation. A DBSIM5D sequence of runs 
was made with step effector deflection angles of 5°, 10°, 

15°, 20 °, 25°, and 30°. The runs were initiated with the 
craft at 56 knts and traveling straight ahead just prior 
to effector deflection. The first peak overshoot value 
and the steady state value attained by the variables 
u (surge) , v (sway) , r (yaw rate) , <j> (roll) , and 0 (pitch) 
were used to document error in the RTS5D. The results of 
the comparison are plotted as error versus effector deflection 
angle in Fig. 18, Fig. 19, Fig. 20, Fig. 21, and Fig. 22. 

B. ACCIDENTAL FAILURE ANALYSIS 

The response of the RTS5D when subjected to accidental 
failure mode scenarios was examined. Two primary situations 
involving the requirement for the pilot to react to the 
failure and maneuver the craft back onto its original track 
were analyzed. The first accidental failure situation 
consisted of reversing the starboard thruster T^ (See 

Fig. 3), such that 6^ equalled 180° while the vehicle was 
tracking straight ahead at 56 knts. The second situation 

consisted of again reversing the T^q thruster to 6^^ equal 
180° with the additional requirement that the failure occur 
while the vehicle was established in a steady state full right 
turn (6_=6 0 =6 f .=6 . =-30° ) . To assume that the "worst case" accidental 

/oy ±U 



72 




6 (degrees ) 

A Steady State Error 
q Transient "First Peak" Error 



Figure 18 

Effector Angle vs Surge Error 
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6 (degrees ) 

& Steady State Error 
O Transient "First Peak" Error 



Figure 19 

Effector Angle vs Sway Error 
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A Steady State Error 
O Transient "First Peak" Error 



Figure 20 

Effector Angle vs Yaw Rate Error 
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A Steady State Error 
Q Transient "First Peak" Error 



Figure 21 

Effector Angle vs Roll Error 
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& Steady State Error 
O Transient "First Peak" Error 



Figure 22 

Effector Angle vs Pitch Error 
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failure scenario prevailed throughout the conduct of the 
two failure situations described above, the pilot was not 
allowed to reduce the thrust magnitudes of either the 
reversed effector T^q or the remaining thrusters T^, Tg , 
and Tg. The pilot was required to recognize the failure had 
occurred, respond with direct real time man generated control 
efforts via the "joystick" (See Fig. 5A) and maneuver the 
vehicle back onto the "north road". 

In each case the exercise was started by initializing 
the craft at the steady state conditions for either a 
straight ahead run or a -30° effector angle right turn. 

The appropriate steady state values for each initialization 
used are shown in Table IV. 



TABLE IV. STEADY STATE CRAFT CONDITIONS PRIOR 
TO ACCIDENTAL FAILURE TEST 



VARIABLE 



STRAIGHT AHEAD 



RIGHT TURN 



u 


87. 


ft/sec 


73.99 


ft/sec 


V 


0 . 


ft/sec 


-7.58 


ft/sec 


R 


0 . 


deg/sec 


2.9 


deg/sec 


<J> 


0 . 


deg 


1.2 


deg 


0 


1 . 


deg 


1.19 


deg 




0 . 


deg 


0 . 


deg 
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Termination of the exercise occurred whenever the pilot 
simultaneously satisfied two recovery requirements: 

1) The craft's navigation coordinates must be within 
50 ft of the north road centerline ( +X NAV axis) . 

2) The translation velocity of the vehicle in the 

NAV 

direction must be less than plus or minus 1 ft/sec. 

It was considered that the craft was completely under control 
once the vehicle had fulfilled the above requirements. It 
is instructional to note that the recovery criteria not 
only required that the pilot counter any asymmetric turning 
moments and sway velocities induced by the failure, but 
that he maneuver the vehicle back to its original track. 

The results of a series of accidental thrust reversal tests 
are shown in Fig. 23A and Fig. 23B. The maximum Y„,. r (See 
Fig. 1) translation was 82.2 ft for the straight ahead run 
case and 571.2 ft for the right turn effector reversal case. 
It was observed that in both cases the craft retained the 
capability to maneuver back to the original track. Further, 
the corrective action required in both cases was simply to 
command the three operational effectors to the full left 
(+30°) position and hold that until the approach to the 
original track, at which time the effectors were modulated. 
Consequently the fastest recovery can be expected to come 
from the operator who is quickest to recognize that an 
effector reversal has occurred. The initial optimal control 
algorithm in both cases was full left rudder. 



79 



100.0 



y nav 




RUN NUMBER 



Figure 23A 
JAV Translal 
Effector Reversal, Straight Ahead Run 



Maximum Y, T .„ Translation For 
NAV 
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y nav 




RUN NUMBER 



Figure 23B 

Maximum Y^y Translation For 
Effector Reversal, Right Turn Run 
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Repeated runs made with the same operator demonstrated 
a very steep learning curve. Rapid improvement in initial 
scores was demonstrated in every case. 
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VI. CONCLUSIONS 



The ability to implement a real-time, man-controlled 
simulator exhibiting the steady state dynamics of the 3K-SES 
has been established. By requiring the right hand side of 
the simplified equations of motion to be restricted to a 
few significant terms, that portion of the iteration loop 
delay dedicated to computation and the AD/DA process has 
been reduced to an acceptable magnitude of approximately 
35 msecs. The remaining portion of the total iteration time 
of approximately 84 msec was required for the subroutine's 
outputting the "alpha-numeric" presentation of the craft's 
states onto a graphics display console. This low value of 
graphics presentation time resulted from the utilization 
of time multiplexing of the data to be displayed. Using 
100 msecs as a maximum allowable iteration time additional 
terms for the right hand side can be incorporated to generate 
a more accurate simulation. From Section III on rectuangular 
integration error analysis it was concluded that the RTS5D 
loop delay must remain less than 100 msec to remain credible. 
This concurs with those conclusions reached by A.L. Greensite 
in Ref. 9. Once the 100 msec loop delay point has been 
reached further expansion of the right hand side of the equa- 
tions of motion would come from development of a more 
efficient graphics system. It was observed that while the 
graphics console utilization required a high cost in terms 
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of loop delay, the interface of realistic man-generated 
control efforts could be accomplished in no other available 
manner due to limitations on existing facilities. 

Results of operational use of the RTS5D indicate the 
craft perspective and north road algorithm developed in 
Section IV-D and used to generate real time visual cues 
for the pilot to be extremely effective in displaying craft 
heading changes and vehicle displacements in Y . 

The results in Section V-A have shown the RTS5D to be 
within acceptable error limitations of less than 2% on the 
average when compared to the DBSIM5D program in the moderate 
speed range of approximately 56 knts and over the full 
effector deflection range of plus or minus 30 degrees. 

The results in Section V-B have shown that the 3K-SES 
is capable of being controlled satisfactorily by man during 
accidental failures involving both loss of thrust and 
effector reversals. From these results it was shown that 
the optimal control effort during the initial stages of an 
effector reversal failure was simply to apply full opposing 
effector angle on the remaining operational effectors. This 
input was held by the pilot until the craft returned to the 
vicinity of the original track at which time the operational 
effectors deflection angles had to be man-manipulated 
slightly to maintain track. Consequently computer assisted 
maneuvering controls would be somewhat helpful but do not 
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appear to be a hard requirement during effector reversal 
scenarios in the 3K-SES. 

It was also concluded that the RTS5D program can be 
readily expanded to a full 6 degree of freedom simulator 
with the addition of simplified equations representing the 
heave motion characteristics of the 3K-SES. 



85 



VII. RECOMMENDATIONS 



The following recommendations are made for future 
studies utilizing the RTS5D to simulate the 3K-SES : 

1. The current right hand side of the craft's equations 
of motion require expansion to include more terms to 
accurately represent the 3K-SES over a wider operating range. 

2. Incorporate the sixth degree of freedom, heave, 
into the simulation. 

3. Include "sea-state" input whenever more computational 
time becomes available. 

4. Program the propulsion plant dynamics to allow 
constant power as an input in lieu of constant thrust. 

5. Design a dedicated micro-processor or mini-computer 
to assist the operator in maneuvering control with the intent 
of more fully examining the computer assisted controls 
requirement in the 3K-SES. 

6. Incorporate actual 3K-SES control hardware into the 
RTS5D to allow more realistic "hands on" testing of not only 
the control hardware itself but also the operator personnel. 
Currently a project is being designed to test operator 
response under various command and failure conditions. 

Results will be reported in the near future. 
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APPENDIX A 





RTS5D PROGRAM NOMENCLATURE 


All 


A/D trunk 500 voltage 


A21 


A/D trunk 501 voltage 


A22 


Added mass coefficient in yaw moment equation 


A31 


A/D trunk 502 voltage 


A33 


Added mass coefficient in roll moment equation 


A3 4 


Added mass coefficient in pitch moment equation 


A41 


A/D trunk 503 voltage 


A51 


A/D trunk 504 voltage 


A61 


A/D trunk 505 voltage 


A71 


A/D trunk 506 voltage 


A81 


A/D trunk 507 voltage 


A91 


A/D trunk 510 voltage 


AASTOP 


performance index limit 


ABSYD 


l-ol 


ABSYOD 


1*0,1 


AD 


array of A/D lines 


A1COA 


speed line index counter 


AM 


mass 


APITCH 


pitch angle in degrees 


APOTK5 


voltage on pot E of thruster console 


APRINT 


# iterations between print execute commands 


AR 


analog turn rate voltage 
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ARATE 


turn rate in degrees 




ARDOT 


analog r voltage 




AROLL 


roll angle 




AS TOP 


performance index limit 




AU 


analog velocity voltage 




AUDOT 


analog u voltage 




AV 


analog sway velocity voltage 




AVDOT 


analog v voltage 




AX 


X-Y plotter scale factor 




AXD 


X-Y plotter X q coordinate 




AXST 


initial condition X-coordinate 


on sea track 


AY 


X-Y plotter scale factor 




AYD 


X-Y plotter Y q coordinate 




AYST 


initial condition Y-coordinate 


on sea track 


BlOl 


D/A trunk T-432 




Bll 


D/A trunk T-420 




Bill 


D/A trunk T-433 




B21 


D/A trunk T-421 




B31 


D/A trunk T-422 




B41 


D/A trunk T-423 




B51 


D/A trunk T-424 




B55 


D/A trunk T-425 




B61 


D/A trunk T-426 




B71 


D/A trunk T-427 




B81 


D/A trunk T-430 




B91 


D/A trunk T-431 




BETA 


drift angle 
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BOAT 


subroutine to generate perspective of SES 


CDP 


bow real forces lumped coefficient 


CDX 


drag forces lumped coefficient 


CDY 


sway forces lumped coefficient 


CDZP 


sidewall rolling moment lumped coefficient 


DAL 


digital to analog call 


DELT 


iteration loop time 


DELTA 


commanded iteration loop time 


DTIME 


v time 


DTIMPLT 


scaled value for v time history 


DXPLOT 


same as DTIMPLT 


DYPLOT 


scaled value for v 


DYREPET 


restart value for v 


EFLG 


error flag for TEXTO call 


ETIME 


r time 


EYREPET 


restart value for r 


EXPLOT 


scaled r time 


EYPLOT 


scaled r 


FI 


stern buoyancy force 


F2 


bow seal pitch force 


F3 


bow buoyancy force 


FSP 


port sidewall buoyancy force 


FSS 


starboard sidewall buoyancy force 


FTIME 


• 

v time 


FTIMPLT 


scaled v time 


FXPLOT 


same as FTIMPLT 


FYPLOT 


• 

scaled v 
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FYREPET 


restart value for v 




G 


gravitational acceleration 




GRAPHO 


subroutine to project image 


(non alpha-numeric) 


GTIME 


R time 




GTIMPLT 


scaled R time 




GXP LOT 


same as GTIMPLT 




GYPLOT 


• 

scaled R 




GYREPET 


restart R value 




HEAD 


craft heading 




HOLD 


subroutine to freeze display/program 


HTIME 


• 

u time 




HTIMPLT 


scaled u time 




HXPLOT 


same as HTIMPLT 




HYBETA 


analog B 




HYBRID 


analog flag coefficient (0 < 


or 1) 


HYPLOT 


analog y axis sea track coordinate 


HYP SI 


analog heading (radians) 




HYR 


analog turn rate (rad) 




HYRDOT 


analog r 




HY REPET 


analog restart y coordinate 


on sea track 


HYSESX 


analog y position of craft 




HYSESY 


analog x position of craft 




HYU 


analog u 




HYUDOT 


analog u 




HYV 


analog v 




HYVDOT 


analog v 
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HYXO 



HYXODOT 

HYYO 

HYYODOT 

ICO 

ICOA 

ICRAFT 

IER 

IIPR 

IIPRD 

IIPRE 

IIPRF 

IIPRG 

IIPRH 

IIPRO 

IIPRP 

IIPRQ 

IIPRR 

IIPRS 

I IX 

IJ 

IJA 

IJAD 

IJAE 

IJAF 

I JAG 

IJAH 



analog X q 
analog X q 
analog Y q 
analog Y q 

multiplex graphics index 
speed line index 
vehicle designator 
error flag 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
time history limit 
print counter 
craft perspective counter 
craft perspective counter 
time history counter 
time history counter 
time history counter 
time history counter 
time history counter 
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IJAO 


time 


history 


counter 


I JAP 


time 


history 


counter 


IJAQ 


time 


history 


counter 


I JAR 


time 


history 


counter 


IJAS 


time 


history 


counter 


IJD 


time 


history 


counter 


IJE 


time 


history 


counter 


IJF 


time 


history 


counter 


IJG 


time 


history 


counter 


IJH 


time 


history 


counter 


IJO 


time 


history 


counter 


IJP 


time 


history 


counter 


IJQ 


time 


history 


counter 


IJR 


time 


history 


counter 


IJS 


time 


history 


counter 



IL 


graphics 


console 


display 


array 


ILA1 


graphics 


console 


display 


array 


ILA2 


graphics 


console 


display 


array 


ILA3 


graphics 


console 


display 


array 


ILA4 


graphics 


console 


display 


array 


ILA5 


graphics 


console 


display 


array 


ILA6 


graphics 


console 


display 


array 


ILA7 


graphics 


console 


display 


array 


ILA8 


graphics 


console 


display 


array 


ILA9 


graphics 


console 


display 


array 


ILA10 


graphics 


console 


display 


array 
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ILA11 

ILA12 

ILA13 

ILA14 

ILA15 

ILA16 

ILA17 

ILA18 

ILA19 

ILA20 

ILA21 

ILA22 

ILA23 

ILA24 

ILA25 

ILA26 

ILA27 

ILA28 

ILA2 9 

I LAP 

IN 

IPACK 

IPERFT 

IPLOT 

IPR 

IPRD 



graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 


graphics 


console 


display 


array 



failure mode counter 

subroutine to load display arrays 

performance index 

graphics console display array 

time history counter limit 

time history counter limit 
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IPRE 


time history counter 


limit 


IPRF 


time history counter 


limit 


IPRG 


time history counter 


limit 


IPRH 


time history counter 


limit 


IPRO 


time history counter 


limit 


IPRP 


time history counter 


limit 


IPRQ 


time history counter 


limit 


IPRR 


time history counter 


limit 


IPRS 


time history counter 


limit 


I SWA 


display array 




ITEXT 


subroutine to display 


• alpha-numeric data 


IVEIWA 


craft perspective/north road array 


IVIEWAA 


sea track array 




IVIEWB 


sea track axis array 




IVIEWC 


u time history array 




IVIEWD 


v time history array 




IVIEWE 


r time history array 




IVIEWF 


v time history array 




IVIEWG 


r time history array 




IVIEWH 


• 

u time history array 




IVIEWO 


p time history array 




IVIEWP 


q time history array 




IVIEWQ 


p time history array 




IVIEWR 


q time history array 




IVIEWS 


analog sea track array 


IVIEWZ 


time history axis array 
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IX 


x axis moment of inertia 


J 


time history counter 


JD 


time history counter 


JE 


time history counter 


JF 


time history counter 


JG 


time history counter 


JH 


time history counter 


JO 


time history counter 


JP 


time history counter 


JQ 


time history counter 


JR 


time history counter 


JS 


time history counter 


L3 


lever arm for buoyancy force 


LD 


draft of sidewall 


LP 


side thrust component lever arm 


LZ 


lever arm for sway drag force 


N 


clock interrupt count 


NAD 


AD trunk line array 


NDA 


DA trunk line array 


00 


lever arm for thrust side component 


OTIME 


p time 


OTIMPLT 


• 

scaled p time 


OXPLOT 


same as OTIMPLT 


OYPLOT 


scaled p 


OYREPET 


• 

restart p value 


P 


roll rate 


PB 


plenum pressure 
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PDOT 


roll acceleration 


PDOTM 


max roll acceleration 


PHI 


roll angle 


PMAX 


max roll rate 


POTK 


thruster console array 


POTS 


control stick array 


POTSA 


function box array 


PSI 


craft heading 


PTIME 


• 

q time 


PTIMPLT 


• 

scaled q time 


PXPLOT 


same as PTIMPLT 


PYPLOT 


• 

scaled q 


PYREPET 


• 

restart q value 


Q 


pitch rate 


QDOT 


pitch acceleration 


QDOTM 


• 

max q 


QMAX 


max q 


QTIME 


p time 


QTIMPLT 


scaled p time 


QXPLOT 


same as QTIMPLT 


QYPLOT 


scaled p value 


QYREPET 


restart p value 


R 


craft turn rate r 


RDOT 


r 


RDOTM 


max r 



READCLOCK subroutine to sample real time clock 
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RESET 


subroutine to reset analog computer 


RHO 


density of water 


RMAX 


maximum yaw rate of craft 


RPERF1 


J 1 


RPERF2 


J 2 


RPERF3 


J 3 


RTIME 


q time 


RTIMPLT 


scaled q time 


RX 


J X 


RXPLOT 


same as RTIMPLT 


RY 


J Y 


RYPLOT 


scaled q value 


RYREPET 


restart q value 


RZ 


x z 


SI 


Tg lever arm 


S2 


T 1(J lever arm 


S3 


Tg lever arm 


S4 


T^ lever arm 


SESX 


scaled X 

o 


SESY 


scaled Y 

o 


SF1 


scale factor in time history arrays 


SF14 


scale factor in time history arrays 


SLIP 


Beta influence coefficient 


SPDLIN 


speedline lower screen limit 


SPEED 


craft's velocity in knts 


STARTCLOCK 


subroutine to start clock 
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STIME 


analog time 


STIMPLT 


scaled analog time 


SXPLOT 


HYSESX 


SYPLOT 


HYSESY 


SXREPET 


restart value for HYSESX 


SYREPET 


restart value for HYSESY 


T 


total thrust 


T10 


thruster number 4 


T10C 


forward vector component of T10 


T10MAX 


max T10 value 


T10S 


side vector component of T10 


T7 


thruster number 1 


T7C 


forward vector component of T7 


T7MAX 


max T7 value 


T7S 


side vector component of T7 


T8 


thruster number 2 


T8C 


forward vector components of T8 


T8MAX 


max value of T8 


T8S 


side vector component of T8 


T9 


thruster number 3 


T9C 


forward vector component of T9 


T9MAX 


max value of T9 


T9S 


side vector component of T9 


TEXTO 


subroutine to display alpha-numeric data 


TFORW 


summation of forward thrust components 
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THETA 


pitch angle 


TIM 


u time 


TIME 


time 


TIMPLT 


scaled TIM 


TINT 


real clock interrupt frequency 


TITLE 


program name display array 


TITLE 0 


program name display array 


TITLE1 


program name display array 


TMAX 


max T 


TSIDE 


sum of side thrust components 


TYAW 


sum of yaw moment 


U 


forward velocity (surge) 


UDOT 


• 

u 


UDOTM 


max u 


UMAX 


max u 


UMUL 


scale factor in speedline generation 


V 


side velocity (sway) 


VCD 


thruster console pot array 


VDOT 


• 

V 


VDOTM 


max v 


VMAX 


max v 


vs ■ 


total velocity 


W 


heave velocity 


WDOT 


heave acceleration 


WE 


width of bow seal 


WRITECLOCK 


subroutine to assign clock interrupts value to 
variable 
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ww 


lever arm for sway drag forces 


X 


Y 

o 


xo 


X 

o 


XODOT 


X 

o 


XDX 


boat perspective scale factor 


XGX 


boat perspective scale factor 


XPLOT 


display array 


XREPET 


reset value for SESX 


XST 


scale factor for analog plot 


YO 


y 0 


YODOT 


Y 

o 


YHIGH 


point D vertical position constant 


YPLOT 


display array 


YREPET 


reset value for SESY 


YST 


scale factor for analog plot 


Z 


operational effector angle 


Z10 


effector number four angle 


Z2 


upper effector angle limitation 


Z3 


lower effector angle limitation 


Z4 


zero angle upper limit 


Z5 


zero angle lower limit 


Z7 


effector number one angle 


Z7I 


modified effector number one angle 


Z8 


effector number two angle 


Z8I 


modified effector number two angle 


Z9 


effector number three angle 
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Z9I 

ZAB 

ZI 

ZZZZ 



modified effector number three angle 
scale factor for Z10 
modified Z 

fixed step effector angle input 
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